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SUMMARY 


Numerical optimization was used in conjunction with an inviscid, full 
potential equation, transonic flow analysis computer code to design an upper 
surface contour for a conventional airfoil to improve its supercritical 
performance. The modified airfoil was tested in the Lockheed-Georgia Com- 
pressible Flow Wind Tunnel. The majority of the test was done at eleven 
million Reynolds number and a four percent tunnel top and bottom wall porosity. 
A limited amount of testing was done to obtain data at other Reynolds numbers 
and wall porosities. 

The modified airfoil's performance was evaluated by comparison with test 
data for the baseline airfoil and for an airfoil developed by optimization of 
only the baseline airfoil's leading edge. While the leading edge modification 
performed as expected, the upper surface re-design did not produce all of the 
expected performance improvements. Although the drag divergence Mach number 
was increased, the modified airfoil exhibited more drag creep than for the 
baseline section. This larger drag creep is attributable to the early formu- 
lation (at approximately Moo=.68) of a relatively strong leading edge shock 
wave . 

Theoretical solutions computed using a viscous, full potential equation 
transonic airfoil code were compared to experimental data for the baseline air- 
foil and the upper surface modification. These correlations showed that the 
theory predicted the baseline airfoil's aerodynamics fairly well, but failed 
to accurately compute results for the upper surface modification. This fail- 
ure is shown to be attributable to the inability of the theory to properly 
treat the thick trailing edge boundary layer associated with the upper surface 
modi f icat ion. 

Numerical optimization is concluded to offer the means for efficiently 
designing advanced airfoils. However, until a completely reliable viscous air- 
foil analysis technique is developed, optimization can be used with confidence 
only when the character of the viscous flow is not s i gn i f icant 1 y a 1 tered dur i ng 
the optimization process. 
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INTRODUCTION 


Efficient transonic performance continues to be an important design 
requirement for many new aircraft. Specialized airfoils whose contours are 
dependent upon design conditions are needed to achieve the desired transonic 
performance. To design these airfoils rapidly and effectively, aerodynamici sts 
must have available accurate and easy-to-use theoretical design methods. The 
nonlinearity of the partial differential equations that describe transonic 
flows has hampered the development of such theoretical design methods. However, 
advances in computational fluid dynamics, together with the availability of 
large and fast computers, have resulted in the recent availability of a number 
of transonic design techniques. 


Airfoil design methods can be categorized as either inverse or direct 
procedures. Inverse methods involve the specification of a desired pressure 
distribution and the calculation of the corresponding airfoil. The need to 
specify a priori a pressure distribution tViat will result in a physically 
realistic optimized airfoil is a disadvantage of inverse procedures. Inverse 
methods have been formulated either by using hodograph equations or by solving 
the problem in the physical plane (e.g., refs. 1 and 2, respectively). Since 
hodograph formulations are applicable only to shock-free flows, they are of 
limited usefulness in transonic design. Also, considerable user expertise is 
required to employ hodograph design methods. Phys i cal -pi ane solutions suffer 
computational difficulties in the leading edge region which are usually avoid- 
ed by specifying the airfoil geometry near the leading edge. Since proper 
leading edge design is needed to design optimized transonic airfoils, this 
approach limits the usefulness of physical plane inverse solutions in transonic 
design work. 


The above-mentioned difficulties are avoided 
In this design approach, a numerical optimization 
a suitable aerodynamic analysis method to design 
sense optimized for specific flight conditions, 
tour can be determined which minimizes drag at a 
with the pitching-moment coefficient constrained 


in direct design methods, 
algorithm is coupled with 
airfoils that are in some 
For example, an airfoil con- 
specific lift coefficient 
to be within acceptable 


limits. Perhaps the most promising direct design method is under development 
at NASA-Ames by Hicks and his associates. The method is described and ^example 
designs are discussed in references 3 and k. Briefly, the method involves 
coupling a numerical optimization scheme developed by Vanderplaats (ref. 5) 
with proven airfoil analysis methods. The ability to use any theoretical 
analysis method and any numerical optimization algorithm makes the technique 
very versatile. 

Currently, the transonic full potential code developed by Jameson (ref. 6) 
is being used to provide the needed aerodynamic data. Viscous effects, known 
to be important in transonic flow calculations, are neglected in Jameson's 
analysis method. An inviscid aerodynamic module is used because currently 
available viscous transonic airfoil analysis programs are not completely re- 
liable, and they require more computation time than inviscid techniques. 
Experience with applying the procedure to design subsonic airfoils has indi- 
cated that aerodynamic performance predicted by inviscid methods will be 
manifested when a viscous analysis of the resulting airfoil is performed. 

There is, nevertheless, a need to experimentally verify the performance of 
transonic airfoils designed using the inviscid procedure. 

Because of its versatility, Hicks' method should be useful both in tran- 
sonic design projects in which performance improvements in existing aircraft 
are sought by airfoil modification, and in projects in which new airfoils for 
advanced aircraft are required. The first application is being investigated 
in a project being conducted by Lockheed. With the aid of Hicks, Lockheed 
researchers recently applied the procedure to the redesign of the forward 12% 
of the C— 1 4 1 airfoil leading edge. This application was selected because that 
airfoil exhibits a large drag creep which might be reduced by a limited modi- 
fication of the leading edge. This work was successful in that the predicted 
drag creep was reduced to the same level obtained by a trial and correction 
process in which analysis methods were used to evaluate many candidate leading- 
edge modifications, but in a fraction of the time. 

The airfoil leading-edge design using Hicks' method was obtained by 
modifying only the C— 1 4 1 airfoil upper surface. In contrast, the airfoils 
designed using the trial-and-correction procedure involved simultaneous upper 
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and lower surface modifications. Consequently, Hicks' modification may be 
easier to incorporate as an aircraft change. Further, because it has a larger 
leading-edge radius, Hicks' airfoil may provide better low-speed stall charac- 
teristics than either the basic C— 1 -4 1 section or the other proposed airfoil 
modifications. 

The Hicks' modification to the C— 1 4 1 airfoil leading-edge, as well as the 
basic C— 1 4 1 section, have been tested in the Lockheed Compressible Flow Wind- 
Tunnel (CFWT) at transonic speeds under an internal Lockheed research program. 
Airfoil surface pressure distributions and the attendant forces and moments 
were obtained for an extensive range of Mach numbers and lift coefficients. 

These data substantiated the theoretically predicted performance improvements 
resulting from numerical optimization. The wind tunnel test data showed that 
a 1 % improvement in (ML/D) may have resulted from the modification of the upper 
surface of the forward 12% of the airfoil. In addition to producing an effic- 
ient airfoil modification, numerical optimization required about half the comp- 
utation time and resulted in a 25% reduction in engineering hours when compared 
to a conventional trai 1 -and-correct ion approach. 

The purpose of the work reported herein was to determine the applicability 
of numerical optimization in an extensive redesign of a conventional airfoil to 
improve its supercritical performance. The problem selected was the re-contour- 
ing of the entire upper surface of the baseline C— 1 4 1 A wing airfoil section. 

This problem was selected because the availability of the work already done on 
the C— 1 4 1 airfoil upper surface leading-edge modification premitted an efficient 
comparison of the use of numerical optimization for limited and extensive air- 
foil modification. 

In this report, the aerodynamic design of the upper surface modification 
using numerical optimization is discussed, the wind tunnel model design and 
test are described, the redesigned airfoil performance is compared with that 
of the baseline and modified leading-edge airfoils, and the design procedure 
is evaluated. 

For reasons which will become apparent, the upper surface modified airfoil 
will be referred to as C 1 4 1 H7472 . The baseline airfoil and the optimized 
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V 


optimization design variables 
Z coordinate normal to airfoil chord line, cm (in.) 

a geometric angle of airfoil chord line, degrees 

y ratio of specific heats 

t wind tunnel wall porosity, % 

Subscripts: 


t . e. trailing edge 

T transition strip location 

0 zero normal force 

1 tunnel station one chord length downstream of model 

oo denotes freestream conditions 

Abbreviations: 

CFWT Lockheed Compressible Flow Wind Tunnel 

1 , 1 . s. lower surface 

u, u.s. upper surface 

L.E. leading edge 


AIRFOIL DESIGN 


Problem Definition 

The design objective of this study was to minimize the cruise drag of the 
baseline C— 1 4 1 airfoil, and to increase the drag divergence Mach number at 
cruise lift. The section Mach number and lift coefficient corresponding to 
airplane cruise conditions are .72 and .57, respectively. In order to have 
a direct comparison with the leading-edge upper surface mod i f i cat ion of the 
baseline airfoil designed by Lockheed and Ames, the airfoil modification was 
restricted to the upper surface. Geometric constraints imposed on the airfoil 
modification were that the th i ckness-to-chord ratio not be reduced and that 
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the airfoil leading-edge retain at least the same degree of bluntness. The 
latter condition was imposed to avoid potentially poor airfoil stall 
characteristics. 

Numerical Optimization 

The airfoil upper surface contour needed to attain 
was determined using the previously discussed NASA-Ames 
optimization scheme. 

Design point definition . - Since the optimization scheme used an inviscid 
transonic method to provide aerodynamic data required during the optimization 
process, the first step in the design study was to determine an inviscid de- 
sign condition. The underlying hypothesis in this inviscid optimization is 
that improvements made at the inviscid design conditions will be realized in 
a viscous flow. 

V 

The inviscid design condition was determined by first computing a viscous 
solution using the NYU transonic flow analysis routine (ref. 1) for the base- 
line airfoil (ClAl-1) at cruise conditions. Since the evaluation of the 
optimized airfoil would be done using wind tunnel data, the solutions were 
computed at the tunnel Reynolds number. Thus, the following conditions were 
specified in performing the viscous calculations: 

C n = -57 

Moo = .72 

R N = 11 x 10 6 
X/C = .10 

The results of these calculations are shown in figure 1 where they are 
compared with experimental data. The agreement between theory and experiment 
shown here is fair, with the major discrepancy being the slightly more aft 
theoretical shock location. Since recent work by Blackwell, et al., reported 
in reference 7 shows the strong dependence of shock location to wall porosity 


the design objective 
aerodynamic numerical 
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and since the theoretical calcualtions are done with free-air boundary con- 
ditions, no attempt was made to seek better theory-to-exper iment agreement. 

The viscous solution indicated that the airfoil angle of attack to 
achieve the above conditions was approximately two degrees. The assumption 
was made that the lift loss due to viscosity for the optimized airfoil would 
be about the same as for the baseline airfoil. The inviscid design point was 
thus defined to be: 

Moo = .72 

o=2° 

Design objective . - The design objective (i.e., the parameter to be 
minimized at the design point) for this study was wave drag minimization at 
the design point. In an attempt to increase Mqq, and at the same time avoid 
a point optimization which would be reflected in a local bucket in the C d 
variation with Mach number, a secondary design objective was to reduce C c j w at 
.74 Mach number. This second design objective was introduced as a constraint 
in the optimization scheme, and its imposition is discussed in the following 
section. Consequently, the objective function in the minimization scheme was: 

OBJ = C dw (a = 2°, Moo = .72) 

Design constraints . - Four constraints were imposed during the optimiza- 
tion scheme: 

(1) Z upper /C (X/C = .5) > .074 

(2) Z upper /C (X/C = .005) > .01 

(3) C n > .85 

(4) C dw (a = 2°, M*. = .74) > .0020 

The first constraint was imposed to ensure that the optimized airfoil was at 
least as thick as the baseline section. An arbitrary nose bluntness was forced 
by the second constraint. The third constraint was required in an attempt to 
ensure that when viscous effects were taken into account, the airfoil would 
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produce at least the desired cruise lift. The final constraint enforces the 
secondary design objective by requiring Cd w at M® = .7^ be less than 20 counts. 
The 20-count level was selected to be compatible with the drag at Moo = .72 to 
produce a flat Cd w versus Moo curve at C n = .57- 

Design variables . - Proper selection of design variables is imperative if 
the design objective is to be efficiently attained in numerical optimization. 

In the leading-edge modification study a single polynomial representation of 
the forward 12% of the upper surface was used. In that case, the design vari- 
ables (i.e., the parameters perturbed during the optimization scheme) were the 
coefficients of the polynomial and/or the exponents on the polynomial terms. 
This approach proved to be successful because only a few terms were required 
to achieve sufficient design flexibility and hence computat ion t i mes were small. 

However, numerous polynomials would be required to provide an entire upper 
surface parameterization with adequate design flexibility with each polynomial 
maintaing ordinate and at least first derivative continuity at the match points. 
Such an approach would not only be complex, it would also be computationally 
expensive. Therefore, an alternative airfoil parameterization scheme was used 
in this study. The scheme was developed by Ames researchers and it involves 
the use of perturbation shape functions to distort the upper surface contour 
of the baseline airfoil. The shape functions used in this work are shown in 
figure 2 with their defining equations. In this case, the design variables 
are the pre-mul t i pi y i ng coefficients which determine the magnitude of the in- 
dividual shape functions. These pre-mul t i pi y i ng (or participation) coefficients 
are adjusted by the optimization scheme until the design objective is met with- 
out violating the constraints. Thus, twelve (12) geometric design variables 
were used in this study. An additional design variable, the angle of attack 
(a), was tried and found to be unnecessary. 

Airfoil nomenca 1 ture . - Since the upper surface modification was developed 
with the objective of reducing the drag of the basel i ne C-1^1 airfoil at .72 
and Mach numbers, and since perturbation shape functions were developed by 
_Hicks at NASA-Ames, the resulting airfoil will be referred to as Cl4lH7^72. 

For the purposes of this report, the baseline section will be called C141-1. 
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Since the optimized leading edge was the sixth in a series of leading edges 
designed for the baseline section, this airfoil is called C 1 4 1 —6 . 

Optimization results . - The numerical optimization was done by starting 
with the baseline airfoil at a = 2° and computing solutions for Moo = .72 and 
Moo = -Ik for the perturbations of the design variables until the design objec- 
tive and the four constraints were met. The initial and final inviscid pres- 
sure distributions for Moo = 0.72 and Moo = ,~]h for the C141-1 and the C 1 4 1 H7^ 72 
airfoils are shown in figure 3- The amelioration of the inviscid flow field 
and the attendant reduction in wave drag resulting from numerical optimization 
is evident in these data. 

The C141-1 and C1^1H7^72 airfoil geometries are shown in figure k and the 
coordinates of the C 1 4 1 H7^72 airfoil are listed in table I. Evident in this 
figure is the attempt by the optimization code to use aft-camber to achieve 
the design objective. Since the lower surface was fixed, the only way to in- 
corporate aft-camber was to "hump" the upper surface at about 75% chord. 

The aft hump leads to a strong adverse pressure gradient on the upper 
surface near the trailing edge which conceivably could have a catastrophic 
effect on the boundary layer. This possibility was examined by computing the 
viscous flow about both airfoils using the NYU 2-D transonic airfoil program. 
The results of these calculations are shown in figure 5 where inviscid and 
viscous results for the airfoils at a = 2° and Mco = . 72 and Mco = .7^ are 
shown. The separation point predicted by the Nash-Macdonal d (ref. 1) criteria 
used in the NYU program are flagged in the pressure distributions. Separation 
is predicted to occur on both airfoils, with the separation point being further 
forward on the C 1 4 1 H 7472 airfoil. The default procedure of the NYU program for 
treating t ra i 1 i ng-edge separations was used in these calculations. (The pro- 
cedure is described in detail in reference 1.) Using the default procedure, 
the program predicts that both airfoils will have about the same drag, despite 
the apparent shock-free flow associated with the modified airfoil. 

The lift loss due to viscosity is also evident in the results shown in 
figure 5, and it is larger than estimated. Consequently, the airfoils have to 
operate at an angle of attack greater than 2° to achieve a C n = -57. 
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The decision was made at this point to test the C 1 4 1 H7^72 airfoil despite 
the predictions of flow separation and the uncertainty in the actual viscous 
design condition. This decision was made because of the belief that test data 
were needed to provi de gu idance for f utu re appl i cat ions of numer ica 1 optimization. 
The model design and wind tunnel test are discussed in the next section of this 
report . 


EXPERIMENTAL TESTS 


Apparatus and Test Procedures 

Mode 1 . - A two-dimensional model of the C141H7472 airfoil was fabricated 
from 17-4PH stainless steel in accordance with the coordinates given in table 
I. The model has a chord of 17-78 cm (7-00 in.) and a span of 50.80 cm (20.00 
in.). The model completely spans the two-dimensional test section and is 
supported from the side walls by means of tangs. A photograph of the model 
installed in the Lockheed Compressible Flow Wind-Tunnel is shown in figure 6. 

The model was instrumented with fifty-three (53) static-pressure measuring 
orifices: 27 on the upper surface and 26 on the lower surface. The orifices 

were located near the mid-span of the model and were mounted flush and normal 
to the local contour. The measured orifice locations are listed in table II. 

The model contour was checked dimensionally using a template and feeler 
gage. Discrepancies from the design contour were mostly within ±.025 mm (.001 
in.) with some areas reaching deviations of ±.050 mm (.002 in.). The airfoil 
surface was polished to conventional transonic model surface finish of less 
than .4 microns (15 microinches). Spanwise twist and warp of the model were 
found to be nil. 

Test fac i 1 i ty . - The model was tested in the Lockheed CFWT. Reference 8 
contains a complete description of this facility. The tunnel is of the blow- 
down type, exhausting directly to the atmosphere. The air storage capability 
is 368 m 3 (13,000 ft 3 ) at 4. 13 x 10 6 N/m 2 (600 psia). A sleeve-type control 
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valve accurately maintains the settling chamber stagnation pressure at selected 
pressure less than or equal to the 1.72 xIO 6 N/m 2 (250 psia) maximum and at 
mass flow rates less than 1 089 kg/sec. (2400 lb/sec.). 

The test section is 50.8 cm (20.0 in.) wide by 71-2 cm (28.0 in.) high 
by 183 cm (72.0 in.) long and is enclosed in a 3-7 m (12.0 ft.) diameter plenum 
chamber. For the model tested, the tunnel height to model chord ratio is 4.0 
and the tunnel span to model chord ratio is 2.9. Model 1 blockage was 3 percent 
of the test section cross sectional area. The top and bottom walls of the two- 
dimensional test section have variable porosity capability (from 0 to 10%), 
obtained by sliding two parallel plates with .635 cm (.250 in.) diameter holes 
slanted 60 degrees from the vertical. The 2-D test section side walls are not 
porous . 

Wake survey rake . - The fixed wake-survey rake used for section drag meas- 
urements is described in figure 7 and shown installed in the tunnel in figure 
6. The wake rake was mounted at the tunnel centerline one chord length behind 
the airfoil model. The rake has a total of 90 total head measurement tubes 
and four static pressure tubes. The tubes are .15 cm (.06 in.) in diameter. 

The wake rake has been calibrated in the tunnel without a model present. 

Data have been obtained in previous CFWT ai rfoi 1 tests similar to that 
conducted herein with the wake rake installed and removed to determine its 
influence on the flow over the airfoil. These unpublished data indicated 
the wake rake had negligible effects on the normal -force coefficient, the 
pitching-moment coefficient, and the airfoil pressure distribution. 

Instrumentat ion . - Measurements of the static pressures on the airfoil 
surfaces and the wake rake pressure were made using electronically actuated 
pressure scanning valves. The full-scale range of the quarter-percent accura- 
cy for the wind tunnel conditions tested: wake rake — ±8.6 dynes/cm 2 (±12.5 

psi); and airfoil pressures ±34.4 dynes/cm 2 (±50.0 psi). CEC force balance 
pressure transducers were used in conjunction with CEC servo amplifiers to pro- 
vide a precise measurement of the atmospheric pressure, stagnation pressure, 
and test section static pressure to .025% of the transducer capability: 6.89 x 

10 s N/m 2 (100.0 psi) for the test section static and 1.38x10 6 N/m 2 (200.0 psi) 
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for stagnation pressure. These transducers allow determination of the test 
section Mach number to a accuracy of ±.001 at the higherst stagnation pressure. 

Angle of attack was measured with a calibrated potentiometer operated by 
the angle of attack drive mechanism. 

Raw pressure data were recorded on magnetic tape utilizing the CFWT high 
speed data acquisition system. 


Test and Methods 

Test conditions . - The aerodynamic characteristics of the C 1 4 1 H 7-4 72 air- 
foil were investigated over a wide range of test conditions. The angle of 
attack of the airfoil was varied from 0 to 5 degrees and the Mach number range 
investigated was from 0.45 to 0.78. Tests were conducted at nominal Reynolds 
numbers based on airfoil chord of 4 and 11 million. The majority of the tests 
was conducted for a wind-tunnel porosity of h%. Limited data were obtained 
over a wall porosity range of 2 to 6% to ascertain the sensitivity of the new 
airfoil to this parameter. These results are included in Appendix A. 


Data reduction . - The static pressure measurements at the airfoil surface 
were reduced to standard pressure coefficients and then machine integrated to 
obtain section normal-force and section pitching-moment coefficients about the 
quarter chord using the following equations: 


Cn 


.1 

(c Pji - Cp u ) d(X/C) 


Cm 


[ (CP* " C Pu) (0 ' 25 ' X/C) d ( x / c ) 

0 


Section 

measurements 


profile drag measurements were computed from the wake survey rake 
by the method of reference 9 utilizing the following equations: 
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The ACj is a correction for the wake rake total head tube displacement 
effect when in a transverse velocity gradient. This correction is discussed 
in reference 9 and is given by 


4C d - 0.36 £ C’ dm . 

Trans i t ion . - For the Reynolds number tests of k million, boundary layer 
transition was fixed. The airfoil was investigated with roughness particles 
located on both surfaces at 0.05 C. The roughness height was 0.00039 C and 
was selected according to the criteria of reference 10. The roughness strips 
were 0.13 cm (0.05 in.) wide and consisted of Ballotini glass beads set in a 
plastic adhesive. 

For the Reynolds number tests of 11 million no transition strips were 
ut i 1 i zed . 

Tunnel wall effects . - The effect of the tunnel walls on the two dimen- 
sionality of the flow is considered to be small. This statement is supported 
by measurements reported in reference 11 on the baseline C-l4l airfoil with a 
similar test arrangement. These results indicated very little variation of 
the pressure coefficient across the model span for various flow conditions. 

The conclusion was further substantiated by observation of oil flow patterns 
at the airfoil-wall intersction. Disturbances in this juncture were confined 
to a very small regions. 

Standard subcritical wind-tunnel boundary corrections (normal -force inter- 
ference and blockage ) have been calculated for this test using the method of 
reference 12. Recent studies (ref. 7) have shown, however, this method to be 
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generally inadequate at transonic speeds. As a result, these corrections have 
not be applied to the data presented herein. 


RESULTS AND DISCUSSION 


A complete set of basic aerodynamic force and moment data as well as 
surface pressure distributions for the C141H7472 airfoil at a Reynolds number 
of 11 million and wind-tunnel wall porosity of b°/ 0 is contained in Appendix A. 
These data will form the basis for assessing the modified airfoil performance 
to be discussed below, and unless stated otherwise all data will be for these 
conditions. Appendix A also contains limited data for a Reynolds number of 
4 million, wind-tunnel wall porosity data from 2 to 6%. 

The evaluation of the performance of the C 1 4 1 H 7^+72 airfoil will be made 
by comparing test data for this airfoil with data obtained in a recent test 
of the C141-1 and C141-6 airfoils. Following this airfoil performance evalu- 
ation, the efficacy of the design procedure will be examined. 


C141H7472 Airfoil Aerodynamics 

The experimental data obtained from the C141H7472 airfoil model in the 
Lockheed-Georg ia CFWT are discussed in this section. This discussion is fol- 
lowed by a review of the aerodynamic forces resulting from the surface 
pressures. 

Airfoil pressures . - The variation with free stream Mach number of the 
surface pressure distributions at the design normal force coefficient of 0.57 
can be seen in the data shown in figure 8. A number of conclusions can be 
drawn from these results. First, the airfoil does not have the aft-loading 
associated with the cusped region of a modern supercritical airfoil. In fact, 
there is a small negative lift region covering the last ten percent chord. 

This lack of aft loading is attributable to the modification of only the upper 
surface. Since the lower surface could not deform, the amount of aft camber 
that could be designed into the airfoil was limited by geometric constraints. 
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The only way to induce any aft-loading was thus to "hump" the rear upper surface; 
a design change which was done by the optimization design code. 


The airfoil has a substantial suction peak which becomes supersonic at 
about 0.6 Mach number. The supersonic region increases in size with increasing 
Mach number, and is terminated by a fairly strong shock wave for Mach numbers 
greater than about 0.64. The maximum shock strength prior to drag divergence 
is reached at a free stream Mach number of 0.68. In this case, the local Mach 
number ahead of the shock is 1.26, and the shock induces a localized separation 
bubble at the foot of the shock. 

The variations with free stream Mach number of shock strength and location 
are summarized in figure 9. The shock movement is orderly, with the shock 
first forming at Moo = 0.60 at about 2 % chord, and moving aft to approximately 
55 % chord at Moo = 0.78. The rapid increase in shock strength at about 0.68 
Mach is evident in this figure. This early formation of a strong shock wave 
can be expected to have detrimental effect on airfoil drag characteristics. 

Airfoil forces . - The airfoil drag coefficient variation with Mach number 
for Cp = 0.57 is shown in figure 10. These data exhibit a substantial drag 
increase in range of Moo = 0.60 to 0.72. This rapid and undesirable drag in- 
crease is no doubt caused by the formation of the relatively strong leading 
edge shock discussed in the preceeding section. From Moo = 0-72 to Ma> = 0.76 
the airfoil drag remains fairly constant. This result is interesting since 
a "flat" Cd vs Mx> curve in this Mach range was one of the airfoil design goals. 
For Mach numbers greater than 0.76, the airfoil experiences the onset of rapid 
drag rise associated with the increasing strength of the shock wave. 

Comparison of Airfoils 

Airfoil C141H7472 was designed at an inviscid condition which was expected 
to produce a minimum wave drag airfoil at a normal-force coefficient of 0.57. 

In this section, the performance of this airfoil at its design point relative 
to both the C141-1 and C141-6 airfoils is evaluated. The variation of meas- 
ured drag with Mach number for the three airfoils at C n = 0.57 is shown in 
figure 11. The following observations can be made form these data: 


1 
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1. Airfoil C141H7472 has a higher drag divergence Mach number than 
either of the other two sections. 

2. For Mach numbers less than Mqq for the C141-1 airfoil, the C141H7472 
airfoil has substantially more drag. 

3. The leading-edge modification to the baseline airfoil (C141-6) reduces 
the supercritical drag and increases Mqq relative to the baseline airfoil 
(C141-1) . 

The reasons for these results can be deduced from the chordwise pressure 
distributions shown in figure 12. For M^ = 0.55, the flows are subcritical, 
and in the absence of any indications of flow separation, the higher drag of 
the modified airfoil is probably attributable to its blunt trailing edge. The 
rapid drag increase of the airfoil C141H7472 at approximately Mo = 0.68 can be 
seen to be due to the formation of a strong leading-edge shock wave which does 
not appear on either of the other two airfoils. The Mach number upstream of 
the shock wave is approximate! y^ 1 . 26 at Mo = 0.68 and C n = 0.56. 

As Mo increases and the design Mach numbers are approached, the pressure 
distributions show that the shock on airfoil C141H7472 moves downstream, and 
decreases in strength. At the same time, a shock wave forms on the other two 
airfoils, and increases rapidly in strength. This behavior results in the 
increased Mgg associated with the C141H7472 airfoil. 

Not directly an outcome of this study, but nevertheless of interest, is 
the performance of the airfoil with the leading-edge modification (C141-6) 
near its design point. When compared with the baseline airfoil (C141-1), the 
leading-edge modification produces significantly more leading-edge suction 
(fig. 12). Also, the suction peak is followed by a nearly isotropic re- 
compression. This behavior results in the avoidance of the premature shock 
formation which occurred when the entire upper surface was modified. An inter- 
esting question, and one which remains to be definitely answered, is the fail- 
ure of the optimization scheme to find a leading-edge geometry which would 
produce a similar isentropic compression when the entire upper surface was 
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modified. Perhaps the different results are attributable to the use of differ- 
ing shape functions in the optimization scheme. 


The airfoil C 7 4 1 H 7^*72 did not perform as expected at C n = 0.57. The 
preliminary theoretical analysis discussed in the airfoil design section in- 
dicated that the normal-force loss due to viscosity might have been under- 
estimated in establishing the inviscid design point. Consequently, the airfoil 
C 1 ^ 1 H 7^+72 might be expected to perform better at a reduced normal-force 
coefficient. However, the drag variations with Mach number for the three air- 
foils at C n = 0-50 exhibit the same general characteristics as they did at 
C n = 0.57, as evidenced by the data on figure 13. 


Design Method Evaluation 

The experimental data comparisons discussed above show that the upper 
surface modification increased Moo at the expense of more drag creep, and that 
in fact the C11*1H7l*72 airfoil's drag at the Mm = 0.72 and Moo = 0.71* design 
points was larger than that of the baseline airfoil. This less than satis- 
factory airfoil performance was not a result of the failure of numerical 
optimization, because the optimization did reduce the design objective, in- 
viscid Cj w . This fact is evidenced by figure 3 where the inviscid pressure 
distributions for both airfoils at the design conditions were shown. 

The lack of performance by the C 1 1* 7 H 7^* 72 airfoil, then, must be due to 
adverse viscous effects which were not taken into account. The effects of 
viscosity on airfoil performance were exami ned by computing solutions using 
the viscous NYU program (ref. 1) for both the baseline and the modified air- 
foils and comparing the solutions with experimental data. The solutions were 
computed using the non-conservative option in the NYU code, and employing the 
program in more-or-less "cook book" form. However, the boundary layer dis- 
placement effects had to be drastically under-relaxed, and two inviscid 
iterations were done between boundary layer calculations to permit a more 
stable convergence process. 

The agreement between theoretical drag predictions and experimental re- 
sults is in general poor for both the baseline and the modified airfoil as 
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demonstrated by the data shown in figure 14. The following observations can 
be made from these data: 

1. Mod ' s predicted to within .00k Mach number for both airfoils. 

2. The shape of the drag curve for the baseline airfoil is fairly well pre- 
dicted, but is underestimated by approximately 10 counts. 

3. Both the level and the shape of the Cl4lH7^72 airfoil's drag curve are 
mispredicted by theory. In particular, the early drag rise occurring at 
M<» = 0.68 is missed in the theoretical calculations. 

The reasons for the failure of the NYU viscous transonic code can be 
deduced by comparing theoretical and experimental pressure distributions at 
the same Mach number and lift coefficient. These data are shown for the 
C 1 1 — 1 and ClAlH7^72 airfoils in figures 15 and 1 6 , respectively. 

The results shown in figure 15 indicate that the C 1 4 1 - 1 airfoil pressure 
distributions are fairly well predicted by theory, with the major discrepancies 
being the shock location, trailing edge pressure recovery, and the lower sur- 
face pressure level. The failure to properly compute the tra i 1 ing edge pressure 
recovery is probably due to the large (approximately 20 degrees) trai 1 ing-edge 
included angle which produces a thick boundary layer. The NYU code uses a 
conventional (albeit adjusted) boundary layer method which is not applicable 
to thick boundary layers. The improper calculation of the trailing edge flow 
causes erroneous results elsewhere, as manifested by the incorrect shock wave 
and lower surface pressures. 

Dispite the fact that the theoretical shock wave is stronger than the 
experimental shock, the NYU program predicts a lower drag level than recorded 
in the test. Drag prediction using this code, then, is probably too uncertain 
for use in numerical optimization, even when the general character of the air- 
foil pressure distribution is fairly well predicted. 

The theoretical pressures for the C 1 4 1 H74 72 airfoil bear little resembl- 
ance to their experimental counterparts, as evidenced by the results shown in 
figure 16. At M® = 0.72, the entire character of the leading edge shock is 
missed, while at Mco = 0.7^, a dual shock is predicted when only one shock 
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occurred. In both cases, the trailing edge pressure recovery is over-predicted; 
even more so than for the baseline airfoil. Also, the lower surface pressures 
are under-predicted, but not quite to the same degree which they were missed 
for the C141-1 airfoil. 

Solutions for the C141H7472 airfoil were also computed at off-design con- 
ditions, and results for Moo = 0.45 and M» = 0.78 are shown in figures 17 and 18, 
respectively (both solutions are for C n = 0.57). For the subcritical Mod = 0.45 
condition, agreement between theory and experiment is good, except for the 
aforementioned trailing edge pressure recovery. At Moo = 0-78, the shock wave 
has moved aft to approximately 55% chord, and the supersonic re-expansion which 
occurred in the Moo = 0.72 and 0.74 solutions is not predicted. Although agree- 
ment here is better than at the design conditions, it is not sufficiently 
accurate for airfoil optimization. In fact, drag predictions are relatively 
accurate only for subcritical flows for the C 1 4 1 H7472 airfoil. 

All of the discrepancies are probably attributable in the main to improper 
modeling of the viscous trailing edge flow. This failure is accentuated for 
the C 1 4 1 H7472 airfoil because the "hump" in the upper surface produces a strong 
adverse pressure gradient near the trailing edge. The resulting gradient pro- 
duces a thicker boundary layer than can be predicted by simple boundary layer 
theory. This thick boundary layer produces a reduced trailing edge pressure 
recovery when compared to the baseline airfoil; a result which is not predict- 
ed by theory. 

Experiences both at Lockheed and elsewhere have shown the NYU code to 
yield reliable transonic results for other airfoils. Consequently, alternate 
reasons for the failure to predict the C141H7472 airfoil's aerodynamics were 
explored. 

One possibility examined was wind tunnel wall interference. Previous 
tests in the CFWT (including the baseline airfoil tests) indicate that 4% 
porosity best simulates free-air conditions. Was it possible, however, that 
good agreement could be attained between theory and experiment for the C 1 4 1 H72 74 
airfoil using a different wall porosity? This question was answered by using 
the limited variable porosity data taken in this test, and comparing them with 
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theoretical solutions at the same lift coefficients. The result of this 
side-study was that good agreement could not be found for any porosity (2% to 
G%) investigated. 

A second possibility for the disagreement might be the failure of the 
inviscid flow region noriconservat i ve formulation. This possibility was ex- 
plored by computing some solutions using the quas i -conservat i ve differencing 
scheme option included in the NYU code. This study was done because a con- 
servative problem formulation is known to be correct for flows with shock 
waves. Consequently there was a need to determine if the use of such a formu- 
lation would improve the data correlations for the modified airfoil. The re- 
sults of this study are summarized in figure 19 where nonconservative and 
quas i -conservat i ve solutions are compared with experimental data for the mod- 
ified airfoil at the design points. As these comparisons show, the theoretical 
solutions are similar, and neither agrees well with test data. The fundamental 
problem of the failure to accurately compute the trailing edge flow is evident 
in both solution techniques. 

An attempt was made to force a match between theory and experiment by 
varying the theoretical angle of attack and/or free stream Mach number. The 
comparisons were made for experimental C n = 0.57 and Me = 0.72, 0.7^*, and O. 78 . 
Theoretical solutions were computed using both conservative and nonconservative 
differencing. The results of that investigation are summarized in figure 20. 
The improved correlation is evident in these data; however, the improvement 
was obtained for substantially different theoretical lift coefficients and 
free stream Mach numbers. Of interest is the close agreement between theoreti- 
cal and experimental angle of attack. This result implies that the good corre- 
lation on shock strength and position was probably obtained because the airfoil 
was at the same effective attitude in both the theoretical calculations and 
experiment. The trailing edge flow is still not correctly predicted, and the 
erroneous drag predictions remain. 

This attainment of the better theory/exper i ment correlation required, of 
course, a priori knowledge of the experimental results. Such a requirement 
obviously cannot be placed on an airfoil design methodology. Hence, the ob- 
servation can be made that an improved transonic viscous airfoil method is 
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required to make airfoil design practical and reliable. The most needed im- 
provement seems to be a better trailing edge flow formulation. 


CONCLUDING COMMENTS 


An upper surface for a conventional airfoil section has been designed 
using numerical optimization to improve the airfoil's supercritical perform- 
ance. The modified airfoil (Cl4lH7^72) was tested in the Lockheed-Georg ia 
Compressible Flow Wind Tunnel. The performance of the modified airfoil was 
evaluated by comparisons with test data for both the baseline conventional 
section (C141-1) and an airfoil development by modification of only the base- 
line airfoil's leading-edge (Cl^tl— 6) . An evaluation of the design procedure 
was then made by comparing theoretically predicted airfoil aerodynamics with 
experimental results. The salient results of the study are summarized below: 

1. The airfoil with modified upper surface (ClAlH7^72) increased Mpp relative 
to the baseline airfoil (C141-1) at the expense of larger drag creep which is 
attributable to the premature formation of a relatively strong shock wave. 

2. Numerical optimization did produce an airfoil with reduced inviscid wave 
drag. However, viscous analysis failed to predict either the premature shock 
formation or the airfoil's drag level. 

3. The failure of the viscous airfoil analysis method in this application is 
probably due to the inability of the method to compute the thick boundary 
layer resulting from the strong adverse pressure gradient occurring over the 
tra i 1 i ng-edge region of airfoil C 1 4 1 H 7^ 72 . 

4. The use of a quas i -conservat i ve formulation in lieu of the standard non- 
conservative scheme did not have a significant effect on the useability of the 
theoretical results. 


5. The concept of numerical optimization offers an efficient and versatile 
method for aerodynamic design. However, inviscid optimization should be 
restricted to limited modifications which do not significantly affect the 
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viscous flow (e.g., the leading-edge re-design study briefly discussed herein), 
or to airfoil designs for which viscous effects are wel 1 -understood* 

6. Research should be devoted to developing an improved two-dimensional 
transonic viscous flow method. If such a method were available, its use 
with numerical optimization would provide a means for the efficient design 
of advanced transonic airfoils. 
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TABLE 1. - DESIGN ORDINATES FOR C141H7472 AIRFOIL 


UPPER SURFACE LOWER SURFACE 


x/c 

Z/C 

X/C 

Z/C 


X/C 

Z/C 

.OOOOO 

.00000 

•35000 

.07070 


.00000 

.00000 

.00020 

.00489 

■37500 

.07185 


.00241 

.00781 

. 00040 

.00625 

.40000 

.07277 


.00961 

.01527 

.00060 

.00724 

.42500 

.07347 


.02153 

.02192 

.00080 

.00804 

.45000 

-07395 


.03806 

-02735 

.00100 

.00874 

.47500 

.07421 


.05904 

.03228 

.00200 

.01133 

.50000 

.07426 


.08427 

.03664 

.00300 

.01320 

.52500 

.07410 


.11349 

.04016 

.00400 

.01471 

.55000 

•07374 


.14645 

.04287 

.01000 

.02069 

•57500 

.07318 


.18280 

.04516 

.02000 

.02653 

.60000 

•07239 


.22221 

.04697 

.03000 

.03059 

.62500 

■07135 


.26430 

.04820 

. 04000 

.03385 

.65000 

.07002 


.30866 

.04895 

.05000 

.03664 

.67500 

.06834 


.35486 

.04928 

.06000 

.03911 

■70000 

.06624 


.40245 

.04891 

.07000 

.04133 

.72500 

.06366 


.45099 

.04784 

.08000 

.04337 

.75000 

.06051 


.50000 

.04619 

.10000 

.04698 

-77500 

.05672 


.54901 

.04398 

.12500 

.05083 

.80000 

.05222 


.59755 

.04093 

.15000 

.05414 

.82500 

.04698 


.64514 

-03717 

.17500 

.05704 

.85000 

.04108 


.69134 

.03314 

.20000 

.05963 

.87500 

.03476 


•73570 

.02927 

.22500 

.06197 

• 90000 

.02840 


•77779 

.02524 

.25000 

.06410 

.92500 

.02240 


.81720 

.02140 

.27500 

.06604 

.95000 

.01673 


•85355 

.01769 

.30000 

.06779 

.96000 

.01446 


.88651 

.01410 

.32500 

.06934 

.97000 
. 98000 
1 .00000 

.01214 

.00978 

.00500 


.91573 
.94096 
.96194 
•97847 
•99039 
•99759 
1 .00000 

.01072 

.00765 

.00498 

.00278 

.00122 

.00030 

.00000 
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Figure 1. Design point 
pressures for 


leoretical and experimental 


:he baseline airfoil (C1*t1-l) 












TERM 

P(1) 

P(2) 

P (3) 

P(J), J=4, 10 
P(11) 


EQUATION 

10(1 - X)X V ( 12 )/e 20X 
10(1 - X)X V ( 12 Ve l,0X 
/~X(1 -X)/e3X 
sin5( X n j) 


WHERE X IS INTERPRETED AS X/C 



EXPONENT TABLE 

j 

n 

k 

.5757166 

5 

.756^708 

6 

1 . 

7 

1 -356915 

8 

1.9^3358 

9 

3.106283 

10 

6.578813 


1 1 

z/c = y 


V(N) P(N) 


Figure 2. Upper surface modification airfoil parameter izat ion 
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Figure 5- Concluded. 
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9- - Effect of Mach number on airfoil shock characteristics at 
C n = • 57, R N = 11 x10 6 , X/C t = free. 
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SYMBOL 

AIRFOIL 

Cn 

Cd 
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C141-1 
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.0092 

-.0433 

0 

Cl 41-6 

.531 

.0090 

-.0406 

□ 

Cl 417472 

.569 

.0096 

-.0417 



x/c 


Figure 12. Comparison of measured airfoil pressures , Rw = 11 xIO 
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Figure 12. Continued. 
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Figure 12. Continued 
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Figure 15. Concluded. 
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SYMBOL DATA MACH a C n C d 


O EXP. .72 3.0 0.57 .0121 

NON-CONS. .72 3.0 0.69 .0096 

QUASI-CONS. .71 3.0 0.65 . 0080 



X/C 

Figure 20. Airfoil ClAlH7 i t72 theoretical and experimental 
pressure match. 
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SYMBOL DATA 

MACH 

a 

^n 

Cd 

O EXP 

0.74 

2.77 

0.57 

.0120 

NON-CONS. 

0.73 

2.75 

0.66 

.0089 

QUASI -CONS. 

0.725 

2.75 

0.68 

.0092 



x/c 

Figure 20. Continued. 
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SYMBOL 

DATA 

MACH 

a 

c n 

c d 

o 

EXP 

OO 

2.5 

0.58 

.0160 


NON-CONS. 

■ 75 

2.4 

0.65 

.0089 

— 

QUASI -CONS. 

.76 

2.5 

0.60 

.0138 



Figure 20. Concluded. 
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APPENDIX A: PLOTTED TEST DATA 

Results from Lockheed-Georgia Comressible Flow Wind Tunnel Test 029 are 
detailed in plotted form in this appendix. The aerodynamic force coefficients 
for airfoil C141H7472 are shown in figures A1 through A9 for the basic test 
conditions of Rjj = 11x10 6 , t = b%, and free transition. Figures A10 through 
A13 contain the aerodynamic coefficients for R|g = 4x10 6 , t = and transition 

fixed at 5 % chord. 

The chordwise pressure distributions are shown for the basic test con- 
ditions in figures Al4 through A29, and for Rfj = 4 x 10 6 in figures A30 through 
A36. Figures A37 and A38 contain the chordwise pressure distributions for wall 
porosities of 2%, 3%, 4%, 5%, and 6%. 
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Figure A2 • - Airfoil force data forH=0.55, Rn = 1 1 x 10 6 , X/Cy = Free 
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Figure A12 • " Airfoil force data for M= O. 76 , R^ = 4x10 6 , X/Cj=.05 
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Figure A17.- Airfoil Pressure Distribution for M = 0.60, = 

X/C-j. = Free 


11 x 10 
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Figure A18 . - A i r fo i 1 Pressu re Distribution for M = 0 . 64, = 1 1 x 1 0 6 , 

X/Cj = Free 


79 





































































88 


Figure A28.- Airfoil Pressure Distribution for M = 0 . 78 , R^-llxlO 6 
X/Cj = Free 
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Figure A30.- Airfoil Pressure Distribution for M = 0.72, R^-^txlO 6 , 
X/C T = .05 
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Figure A31 Airfoil Pressure Distribution for M = 0.74, •RN = Ztx1 ° 6 . 
X/C T = .05 
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Figure A34 


Airfoil Pressure Distribution for M = O. 76 , R^ - ** x 10 , 

x/c T = .05 
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Figure A36-" Airfoil Pressure Distribution for M = 0.78, R.. 

X/C T ■ .05 " 
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